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Abstract 

The distribution of finite time observable averages and transport in low dimen¬ 
sional Hamiltonian systems is studied. Finite time observable average distribu¬ 
tions are computed, from which an exponent a characteristic of how the max¬ 
imum of the distributions scales with time is extracted. To link this exponent 
to transport properties, the characteristic exponent /r(g) of the time evolution 
of the different moments of order q related to transport are computed. As a 
testbed for our study the standard map is used. The stochasticity parameter K 
is chosen so that either phase space is mixed with a chaotic sea and islands of 
stability or with only a chaotic sea. Our observations lead to a proposition of a 
law relating the slope in 5 = 0 of the function /^(g) with the exponent a. 


1. Introduction 

The question of transport in Hamiltonian systems is a long standing issue as 
it can be inferred from the vast literature on the matter and references therein [ll 
nil [sum [ 7 ] and the large domain of applications ranging from hot magnetized 
5 plasmas to astronomy, chaotic advection, underwater acoustics etc... Beyond the 
fully chaotic situation in which we usually can apply the central limit theorem, 
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and therefore still have a random walker picture in mind, problems are still not 
clear when the phase space is mixed. 

Indeed in this situation the system is not ergodic, in the sense that there 
10 is not only one unique ergodic component, but instead there are regions with 
chaos, and regions with regular motion. When considering one-and half degrees 
of freedom system, one usually talks about a picture with a stochastic sea and 
islands/regions of regular motion. This coexistence can lead to some problems 
especially since it is possible for Hamiltonian systems to have so called sticky 
15 islands. This paper inscribes itself in this series and tries to tackle the prob¬ 
lem of transport using distributions of finite time observable averages and their 
evolution as the average is computed over larger and larger times. A first at¬ 
tempt using this approach was performed in [ 5 ], using a perturbed pendulum 
as a case study, and it found its roots in the study of advection described in [ 9 ] , 
20 where finite time averages were used to detect sticky parts of trajectories. In 
the present case we simply use the standard or Taylor-Chirikov map [T]. This 
choice was motivated by the fact that such maps can be directly computed from 
the flow of the so-called kicked rotor, and, being a map, it allows us to perform 
fast numerical simulations and gather enough data to have somewhat reliable 
25 statistics. The purpose of this paper is not a thorough study of transport in the 
standard map, but to use this map as a testbed for our analysis of finite time 
observable averages. 

Regarding the problem of transport in Hamiltonian systems, the standard 
map has become over the years a classical case study. One of its advantages is 
30 that it depends on just one control parameter K and many attempts were made 
to find the link between K and a diffusion coefficient [n muni ED. Depending 
on the values of K, we can get a system which is very close to an integrable one 
or one that is fully chaotic, with, in between, the picture of a mixed phase space 
with a chaotic sea and regular islands. In this last setting, we can face the so 
35 called stickiness phenomenon: particle’s trajectory originating from the chaotic 
sea can stay (stick) for arbitrary large times in the vicinity of a stable region. 
This type of phenomenon is able to generate long memory effects, which, in turn. 
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can generate so called anomalous transport also called anomalous diffusion. In 
contrast with normal diffusion, the dynamics leads to transport properties which 
40 can be far from the Gaussian-like processes, and the second moment grows 
nonlinearly in time. 

In the following section, we briefly introduce the standard map and present 
the phase space and first results with the choice of parameters we considered. 
Then we discuss and consider transport properties in each system. We present 
45 the method and compute characteristic transport exponents. We confirm the 
multi-fractal nature of transport in both two considered cases where the phase 
space is mixed, while transport appears as diffusive in the global chaotic case. 
Finally, we investigate the relation between a, the characteristic exponent of the 
evolution of the maximum of the distribution of finite-time observable averages, 
50 and /i, the characteristic exponent of the second moment of transport associated 
to the observable. In [5], a simple law was proposed, namely a = 1 — /r/2, 
our findings lead to good agreements for two out of the three cases. As a 
consequence, a slightly more general law is then proposed which captures all 
features; and then we conclude. 

55 2. The standard map 

Before moving to more details, we remind the reader that the standard 
map arises naturally as a Poincare mapping of the kicked rotor model, whose 
Hamiltonian writes 

2 

i? = y - cos(g) X S{t-nT), (1) 

n— — oo 

where the parameters wq and r are without dimensionality and <5 is the Dirac 
function. We shall not derive the standard map here, and we will consider it on 
the torus. In this case its equations are 

{ p„+i =p„ + K sin q„ [27r] 

9n+l = Qn+Pn+l pTf] , 
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Figure 1: Phase space visualization for k = 1.5. We observe a mixed phase space and Hamil¬ 
tonian chaos which covers about half of the phase space. 



Figure 2: Phase space visualization for k = 2.5. We observe a mixed phase space and Hamil¬ 
tonian chaos which covers a larger portion than in Fig. ^ 
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Figure 3: Phase space visualization for k = 10. We observe a fully chaotic phase space and 
no regular islands subsist. 


where K is the parameter that characterizes the force amplitude [T]. Before 
going to study transport in this system, we shall briefly present the three differ¬ 
ent cases considered. Namely, we considered three different values for K. The 
60 plots for K = 1.5, K = 2.5, and = 10 are represented respectively in Figs. 

|1|2[ and[^ 

Now that we specified the object of our study, let us consider transport. 

3. Transport Properties 

In order to consider transport, we shall first consider an observable. In 
previous studies [HI 11 [la ng, it has been found that considering the absolute 
speed as an observable, leads to relatively clear results for transport. However 
in the previously mentioned works, the considered Hamiltonian systems were 
flows, and the study of transport was made by considering the dispersion of 
the arc-length of different trajectories. Here we consider the standard map 
<§■ The notion of arc-length does not make much sense unless we consider the 
underlying kicked-rotor flow. Though, we may consider an analog to the norm 
of the phase space speed, by considering the distance between the two points 
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(9n+l tPu+i) and {qn,Pn)', given the equations (|^. This unfortunately does not 
define a proper observable since it mixes both steps n and n + 1. In order to 
get a real phase space observable we define one inspired from this distance by 

V = \J K'^ sin^ q + p"^ . (3) 

65 We shall now consider the average of v along a typical trajectory. If the 
system is ergodic we may naturally expect that the average will converge to 
the phase space average of v with the ergodic measure. In order to assess 
this statement we shall consider initial conditions in the stochastic sea, and 
follow trajectories for large times; and following the ideas developed in mm, 
70 we consider the distribution of finite-time averages. This means that we shall 
compute averages of v over finite times, namely for the given initial condition 
(9i, Pi) we compute 

^ n—1 

Vzin) = , (4) 

k=0 

where Vi(k) is the value of the observable v at the step k for a trajectory given 
by its initial condition pi). Let us denote Pniv) the probability density of the 
distribution of Vi{n)'s Q. Assuming ergodicity and we introduce the ergodic 
average (w). We then have 


lim p„(v) = 5(1; - (n)) . (5) 

n—^oo 

Since the Dirac function is singular, one simple way to asses the convergence of 
Pm assuming it is relatively smooth, is to just look at how fast its maximum 
value pmax{n) grows towards oo with n (see for instance Fig. top plot). We 
shall see later that we have to be careful when defining Pmaxin), but it is easy 
to figure out that should the dynamics be sufficiently chaotic, so that a central 
limit theorem applies, we can genuinely expect 

Pmax{n) ~ n“ = Vn , (6) 
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Figure 4: fc = 10. Top: distributions of for from flattest to thinnest n = 10^, 8 x 
10^, 3.2 X 10“^, 1.28 X 10^. Bottom: evolution of pmax{n) versus n, in logarithmic scale. One 
can observe a scaling pmax ~ implying regular diffusive transport with ol = 1/2. Note 

that the average has been removed so that distribution are all centered on zero. Data is 
taken from 1024 trajectories with initial conditions in the stochastic sea and computed for 
10® iterations. 
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Figure 5: fc = 2.5. Top: distributions of Un, for from flattest to thinnest n = 10^, 8x 10^, 3.2 x 
10*^, 1.28x10^. Bottom: evolution of p 7 naa;(^) versus n, in logarithmic scale. One can observe 
a scaling pmax ~ with a ^ 0.43 superdiffusive transport. Note that the average has been 
removed so that distribution are all centered on zero. Data is taken from 1024 trajectories 
with initial conditions in the stochastic sea and computed for 10^ iterations. 


































r 




' 












_ 


0 1 - 1 - 1 - 1 -^- 1 

3 4 5 6 7 8 

logi()(^) 


Figure 6: fc = 1.5. Top: distributions of Vn, from flattest to thinnest distribution n = 
10^, 8 X 10^, 3.2 X 10^, 1.28 x 10^. Bottom: evolution of pmax{n) versus n, in logarithmic 
scale. One can observe a scaling pmax ~ with a ^ 0.38 implying as well superdiffusive 
transport. Note that the average has been removed so that distribution are all centered on 
zero. Also there are strong thin peaks due to stickiness in the distributions, and that pmax 
has been computed by measuring the local maximum of the central flat peak located near 
Vn = 0. Data is taken from 1024 trajectories with initial conditions in the stochastic sea and 
computed for 10® iterations. 
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where we defined an exponent a which is equal to 1/2 in this case. To convince 
ourselves we may consider the variables 


rj — 1 



whose distribution converges towards a Gaussian when the central limit theorem 
applies. Since between j/’s and v’s, there is ’’just” a rescaling by \l\fn, one 
expects that the mean square displacement of v shrinks as Then, since 

the total area is conserved and equal to one [p is a probability density function), 
it is natural to infer that pn looks more and more like a Gaussian with a variance 
that shrinks as \/n and a maximum that grows therefore as Pmax{n-) ~ ^/n. 
In fact the convergence towards the Gaussian depends on the considered point. 
Here we assume that there are no problems with large deviations and that for 
large enough n everything is under control (see for instance mm)- 

Before moving on to the specific results obtained for the considered cases, 
let us emphasize a last feature. For the sake of analogy we previous studies, we 
shall define what is the arc-length equivalent of the flows for this map as 


n—1 


Siin) = ^ Vi{k) = nvi{n) = ^/nyi{n) -I- n{v) . 
k=0 


(7) 


It is evidently abusive to denote this as an arclength. Anyhow the origin of 
anomalous transport is due to the breaking of the central limit theorem. For our 
case, this breaking can only occur due to strong time-correlations, i.e memory 
effects. Since the central limit theorem in dynamical system is applied to a 
given observable, meaning a function associating a point in phase space to a real 
number we have to define one. In previous works in low dimensional Hamiltonian 
flows a suitable choice for an observable appeared to be the norm of the speed in 
phase space. Using this observable attached to no particular coordinate system 
ended up as a good choice, especially when the phase space was bounded. When 
dealing with the standard map, the notion of speed per se has no meaning, so 
for the sake of analogy with previous work we here choose a specific observable 
V which “resembles” a speed, and its associated displacement s. It is important 


to remind the reader that, except of a few specific observables, the nature of 
95 transport will not depend on the choice of the observable. 

We now consider the transport properties using s. If the central limit the¬ 
orem applies, one will expect that the varianceof the distribution of s^’s grows 
like y/n^ so that the second moment has a characteristic exponent /r = 1. 

Let us now envision a situation for which transport is anomalous, for instance 
super-diffusive, with so called fat tails like power law decreasing ones for instance 
giving rise to a transport exponent 1 < /r < 2. The scaling relation Q between 
s and V still holds, which means that the variance of v decreases not as fast 
as for the Gaussian case, since we still have area conservation, we can expect 
that the maximum grows slower than y/n (if the distribution is flat enough at 
the top) and thus that the characteristic exponent related to the growth of the 
maximum of the distribution a is such that a < 1/2. In fact since, we will 
assume that these power law behaviors of the maximum and the variance are 
valid. Then since the total probability is conserved and equal to one for each 
time n, a fact that can be viewed as area conservation under the drawing of 
the function p„(w). Since we are dealing with scaling laws, we can make the 
rather crude approximation that the probability density function Pniv) can be 
approximated by a rectangle function. Then considering the conservation of the 
area S = pmax{n) x y/ cr(n), where a(n) is the variance of the Vi{n), we directly 
obtain a relation between p, and a (see [5] ) which writes 

a = 1 - I . (8) 

In any case even if Eq.(|^ does not hold, the previous considerations offer a 
100 ’’different” possible way to characterize anomalous transport depending on the 
values of a: 

• If a 1/2 transport is anomalous, moreover 

— if a > 1/2 we expect sub-diffusive transport 

— if a < 1/2 we expect superdiffusive transport. 
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In order to confirm these statements for the considered system ([^, we as well 
compute its transport properties. In order to characterize transport, we follow 
m and consider the transport properties related to the chosen observable, 
meaning we compute the different moments of Q, from which we extract a 
characteristic exponent 

Mq{n) = (|sj(n) - (si(n))|'J) ~ (9) 

105 The second order exponent /r = /i(2) characterizes the anomalous and diffusive 
transport. We recall that transport is so-called diffusive if /i(2) = 1 , and in 
all other cases the transport is anomalous. More specifically, if /i(2) < 1 the 
transport is sub-diffusive and if /r(2) > 1, it is super-diffusive. 

In order to check all these assumptions we considered the standard map 
110 dynamics in the three specific case described in Figs. We numerically 

computed the evolution of 1024 different trajectories, for which we computed 
10® steps. Initial conditions where taken in a square of side O.OI centered on 
the point {p, q) = (0.1, 0.1). Records of s’s and positions are taken every 10® 
steps, meaning that we have about 10® records to compute our distribution and 
115 moments. When computing histograms to obtain the distribution, we sampled 
equally spaced 5 10^ bins between the minimal and maximal value of the p(u)’s. 
We insist on the fact that when binning is not adequate we may end up for large 
times having all the Vi{n) accumulating in the same interval, when this happens 
the distribution we obtain is effectively a rectangle function, whose maximum 
120 does not grow, a phenomenon that can be monitored using Pmax(n). We may 
also end up having problems with the accuracy of our data due to the finite 
precision of our real numbers (we considered here double precision numbers 
with 16 characteristic digits), which would lead to a similar phenomenon but 
may as well accelerate the growth of Pmax{n). Finally if our sample is to refined 
125 in regards to our amount of data, the distribution becomes too discontinuous to 
mean anything and we can not rely on our analysis. Given these constraints, we 
settled for the aforementioned numbers, and checked the stability of the results 
versus different binning strategies. Increasing the duration of the trajectories 


12 






Figure 7: k = 1.5. Top: moments versus time of s. Bottom: characteristic exponent versus 
moment order. These exponent have been computed using the linear portion of the top 
picture. This nonlinear behavior is typical of strong anomalous transport. We have super- 
diffusive transport and fi{2) ^ 1.3. 
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Figure S: k = 2.5. Top moments versus time of s. Bottom characteristic exponent versus 
moment order. These exponent have been computed using the linear portion of the top 
picture. This nonlinear behavior is typical of strong anomalous transport. We have super- 
diffusive transport and fi{2) ^ 1.67. 
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Figure 9: k = 10. Top: moments versus time of s. Bottom: characteristic exponent versus 
moment order. We observe the linear behavior with slope 1/2 characteristic of diffusive 
transport. 
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(computing more time steps) will lead to problems of accuracy due to finite 
130 precision, and it would thus necessitate the use of computing trajectories with 
higher precision. The other strategy would be to compute more trajectories, 
i.e take more initial conditions, however this leads to a bias because it would 
only work assuming that the ergodic measure in the stochastic sea is related 
to the Lebesgue one, and to avoid this bias we prefer to consider finite time 
135 portions of trajectories computed over large time. Of course using more powerful 
computers, computing more data is always possible, but the presented results 
are sufficiently accurate in order to confirm our conclusion on the nature of 
transport and to link between the exponents. 

We now return to the analysis of our results. 

140 The first simpler case is the fully chaotic one, namely the case for which 
K = IQ (Fig. [^. We first represented in Fig. the distributions resulting 
from different time averages, and the evolution of the maximum Pmax versus 
the length over which we average n. We find as anticipated that a = 1/2 
and thus expect to have normal diffusive transport. We then move on to the 
145 two other cases with mixed phase space, meaning with a non-unique ergodic 
measures. For these, we consider initial conditions in the so-called stochastic 
sea. Should we consider trajectories inside the regular regions they would have a 
ballistic like contribution, and would not belong to the same ergodic component. 
For the situation displayed in Fig. we can notice in the distribution that 
150 small peaks (bumps) appear near « 0.3, these peaks correspond to the 
stickiness phenomenon and portions of trajectories that remain a long time 
around a regular island, giving rise to memory effect and Levi flights (see for 
instance m- The presence of these peaks naturally affects a, and we find 
a « 0.43, meaning we expect super-diffusive transport. For the last considered 
155 case K = 1.5 displayed in Fig. we observe a similar phenomenon, the peaks 
are however much sharper than in the K = 2.5 situation. In fact for small 
values of n the presence of these peaks affects the value of Pmax{n), in order 
to circumvent this problem we considered in Fig. only the local maximum 
centered around the final average value; for large values of n this does not affect 
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160 the value of Pmaxin), but it does for small n’s. In this case, the polynomial 
behavior appears to be only roughly correct, and we measure a « 0.38, so we 
expect as well super-diffusive transport. If we have some faith in Eq.([^, we 
expect as well a larger second characteristic exponent for transport than for the 
K = 2.5 case. 

165 In order to verify our conclusions we computed the transport properties from 
the different data sets. Meaning we computed the different moments Mq{n) 
and extracted from each a typical characteristic exponent The results are 

displayed in Figs. [3[8]and|| 

We shall discuss these results starting from the simplest case, namely the 
170 fully chaotic situation corresponding to K = 10. The results of this case is rep¬ 
resented on Fig. as expected in this testbed situation we recover the features 
of Gaussian diffusive transport, meaning a second moment /i(2) = 1, and a lin¬ 
ear behavior of the moments p{q) = q/2 |16| . Given the previous results and 
the value of a = 1/2 measured previously, we have a first setting for which 
175 the equation Fq.([^ holds. Gonsidering a more complicated case, namely when 
K = 1.5 (Fig. [^, the nonlinear behavior of the function /i(q) with /r ~ g for 
large values of q. According to the definitions given in m, this behavior im¬ 
plies that we are facing multi-fractal transport or so called strong anomalous 
transport. Anyhow, we still find a good agreement between the measured value 
180 of the second moment /i(2) « 1.3 and what the value of we computed as well 
a = 0.38, the equation Fq.Q still holds and we insist as well that since a < 1/2 
we found as expected super-diffusive transport. 

It is in fact for the intermediate case K = 2.5, with the small peaks that we 
got unexpected results. Indeed, in Fig.[^ we can notice that the second moment 
185 exponent is located in the area corresponding to the change of asymptote be¬ 
tween two linear regimes. The measured value of the transport exponent gives 
/r(2) « 1.67. This measurement allows to clearly state that the equation Fq. (§ 
is not general enough and does not hold for this particular situation. Indeed we 
measured a « 0.43, while should Fq.([^ be true we would get Oexpected = 0.17. 
190 Note though that the measured value of a correctly concludes on the nature of 
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the transport which is also super-diffusive. 

Measuring this exponent is therefore a good indicator of the nature of trans¬ 
port, but how is it related to the measured transport exponent. In fact in all 
of our previous computations, as well as in the results presented in [5], Eq.(|^ 
195 was holding. In order to explain what is going on we shall have a closer look 
at the exponent versus moment order figures displayed in Figs. mi and[^ We 
can notice that both cases for which the expression Eq.(|^ holds, correspond 
to an exponent whose value can be obtained by simply extrapolating the linear 
behavior of the function ^,{q) for small values of 9 , i.e the small moments linear 
200 behavior. This feature was also true in the results presented in [5], where the 
expression Eq. was first derived. A natural simple conclusion is then to con¬ 
clude that the same is true for this K = 2.5 problematic case. In order to check 
this we perform the same extrapolation from the low values of g to g = 2 of /r(g) 
for the case K = 2.5 in Fig. We observe as well a linear behavior in Fig. |^for 
205 low moments which we can extrapolate, we then measure ^extrapoiated{‘^) ~ 1 - 2 . 
This extrapolated value is in remarkable agreement with what we would expect 
using formula indeed we would then obtain a = 0.4 close to the measured 
a = 0.43 one. 

From these last remarks and measurement we may decide that this linear 
210 extrapolation is actually a more general expression relating the exponent a 
to the transport moments /r(g). We may therefore speculate that a correct 
expression is 




( 10 ) 


meaning that the exponent a is directly related to the slope at the origin of 
the function /r(g). We conjecture, that in all situations for which the previous 
215 expression Eq.® holds, the observed value /r(2) can be extrapolated from the 


slope at the origin. Eq. (10 1 provides therefore a more general setting that is 
compatible with previous results, but also when facing situations like the one 
for K = 2.5. 
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In order to understand why actually this Eq. ( |Io| ) is probably more sound, 
220 we are looking again at the moments displayed in Figs[^and[^ In these pictures 
we notice that the evolution is not of the power law type for large times, with 
some not standard more or less flat evolution. In fact when looking at the data, 
this behavior has a simple explanation. In fact this may be due to the finite 
time sampling and finite number of trajectories. Indeed let us consider the 
225 situation of a trajectory captured in a long Levy flight (captured in a sticky set 
for instance) which suddenly stops ’’flying” (goes back to the stochastic sea). 
When computing high order moments, we can naturally state that these will be 
essentially dominated by this ballistic excursion, and when the flight stops, the 
moments will then remain more or less constant until the rest of the trajectories 
230 finally manage to reach such large excursion. On the other side this effect is 
less visible on lower moments. Indeed the inversion of convexity for q < \ 
puts less emphasis on large excursions. The same is true for the exponent a. 
Indeed the portion of a trajectory caught in a flight will contribute to a specific 
peak different than the bulk one (see for instance Fig. [^, then when the flight 
235 stops, it will take a while before the time average of the speed of this portion of 
trajectory enters the bulk, and starts contributing to the main peak and thus 
influence the growth of its tip, i.e the value of a. 

4 . Conclusion 

In this paper we studied transport properties using distributions of finite 
240 time observable averages and monitoring the time evolution of its maximum in 
the spirit of the work presented in [5]. In order to have access to large amount 
of data, we opted to reduce the numerical computation of the flow, while still 
considering a Hamiltonian system; we opted for the kicked rotor system which 
can be reduced to the standard map, which we settled for as a test case study. 
245 In this setting we have shown that a, the characteristic exponent of the time 
evolution of the maximum of finite time averaged observable distributions, con¬ 
veys good information relative to the nature of transport. Indeed, we confirmed 
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the nature of transport with respect to the values of a, and we found super- 
diffusive transport when a < 1/2 which corresponds to situations of a mixed 
250 phase space, while we recovered Gaussian transport for the fully chaotic regime 
and a = 1/2. We as well confirmed the multi-fractal nature of transport in 
the standard map for the cases with mixed phase space already discussed in 
[Is] and reference therein. Finally we propose a link between the characteristic 
exponent of transport moments and a, which generalizes the results obtained 
255 in [5]. This expression appears to encompass as well situations when the expo¬ 
nent of the second moment lies in the nonlinear zone of the curve relating the 
characteristic moment exponent ^{q) versus the order of the moment q. Using 
this new relation, we find a good agreement between the measured values of a 
and the behavior of the characteristic transport exponent for low order. Finally 
260 one could ask the reason on why to introduce this a exponent. Besides its in¬ 
trinsic interest we found in the end that it is computationally easier and faster 
to compute. Indeed, as long as good care of histogram computation and bin 
size is performed, the linear scaling law to extract its value appears to be much 
easier than when dealing with transport for which this is not always obvious 
265 especially when transport is super-diffusive. 
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